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An overview of the known Wolf-Rayet (WR) population of the Milky Way is presented, 
including a brief overview of historical catalogues and recent advances based on infrared pho¬ 
tometric and spectroscopic observations resulting in the current census of 642 (vl.l3 online 
catalogue). The observed distribution of WR stars is considered with respect to known star 
clusters, given that <20% of WR stars in the disk are located in clusters. WN stars outnumber 
WC stars at all galactocentric radii, while early-type WC stars are strongly biased against the 
inner Milky Way. Finally, recent estimates of the global WR population in the Milky Way are 
reassessed, with 1,200±100 estimated, such that the current census may be 50% complete. A 
characteristic WR lifetime of 0.25 Myr is inferred for an initial mass threshold of 25 Mq. 


1 Historical overview and current 
census 

Wolf-Rayet (WR) stars are the evolved descendents 
of massive stars, with two main flavours: Helium- 
rich WN stars displaying the products of core H 
burning and Carbon-rich WC stars displaying the 
products of core He burning (Crowther, 2007). An 
overview of the number of WR stars in the Milky 
Way is provided. The Galactic census tripled be¬ 
tween the 1st (Campbell, 1894) and the 6th cata¬ 
logues (van der Hucht 1981), a century later (Ta¬ 
ble 1), and subsequently doubled by the Annex 
to the 7th catalogue (van der Hucht, 2006) owing 
to deep narrow-band optical surveys (Shara et ah, 
1999), the discovery of rich star clusters (Crowther 
et ah, 2006) and infrared surveys of the Galactic 
Centre region (Figer et ah, 2002). Over the past 
decade, the number has doubled again, primarily 
through surveys arising from the advent of large 
format infrared detectors, with the census totalling 
642^, comprising 357 WN stars, 8 WN/C stars, 273 
WC stars and 4 WO stars. 


Tab. 1: Historical catalogues of Galactic WR stars 


Catalogue Reference 

Number 

I 

Campbell (1894) 

55 

H 

Fleming & Pickering (1912) 

108 

HI 

Payne (1930) 

92 

IV 

Roberts (1962) 

123 

V 

Smith (1968) 

127 

VI 

van der Hucht et al. (1981) 

157 

VH 

van der Hucht (2001) 

227 

VH 

Current census vl.l3 

642^ 


Two complementary initiatives have dominated 
this improvement since the annex to the 7th cata¬ 
logue (van der Hucht, 2006). Mike Shara and collab¬ 
orators have employed narrow-band imaging of the 


Galactic plane at 2-2.3/im to identify 150+ WR stars 
from follow-up spectroscopy (Shara et ah, 2012). 
Since this approach relies upon the contrast between 
emission lines and the adjacent continuum, systems 
which are dominated by hot dust or a companion in 
the K-band will not be identihed. A second method 
exploiting the unusual near- to mid-IR broad-band 
colours of WR stars has been exploited by Schuyler 
Van Dyk, Pat Morris and collaborators has identi¬ 
hed 100+ WR stars (Hadheld et ah, 2007; Mauerhan 
et ah, 2011). This approach has successfully iden¬ 
tihed dusty WR stars, but would also prove prob¬ 
lematic for systems whose IR energy distribution is 
dominated by a companion star. 

# Optical # Gal Centre # IR narrow-band 

# IR SED # IR clusters # IR nebulae 



Fig. 1 : Pie chart illustrating breakdown of discovery 
technique for WR stars in current census (vl.l3 of on¬ 
line catalogue) 


In addition, WR stars have been discovered 
from follow-up spectroscopy of members in star 


^vl.l3 of online WR catalogue http://pacrowther.staff.shef.ac.uk/WRcat/ 
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clusters newly identified from IR surveys (e.g. 
Spitzer/GLIMPSE, VISTA/VVV) such as Chene 
et al. (2012). Other WR stars have been newly 
identified from follow-up spectroscopy of the ioniz¬ 
ing stars in dusty mid-IR nebulae (Wachter et ah, 
2010). A chart illustrating the fraction of WR stars 
discovered by different techniques is illustrated in 
Fig. 1. 

The binary frequency of O stars in young Galactic 
clusters is believed to be high: ^71% will interact 
with a companion during their lifetime according to 
Sana et al. (2012). Since WR stars are the progeny 
of O stars, a high WR binary frequency will also be 
anticipated. The majority of searches for binarity in 
WR stars have been conducted via blue visual spec¬ 
troscopic monitoring. Amongst the 25% of WR stars 
which are visually brighter than V=15 mag, 44 are 
close binaries (either SBl or SB2), plus another 10 
are dusty WG stars. It is now widely accepted that 
dusty WG stars are universally binary systems in¬ 
volving a WG star and an O star companion. There¬ 
fore, a strict lower limit to the binary frequency is 
34% (= 54/160 stars). For comparison, a WR binary 
frequency of 38% (86/227) was estimated in the 7th 
WR catalogue (van der Hucht, 2001), including long 
period systems (P > 1000 days) and X-ray bright 
systems - by way of example the X-ray bright sys¬ 
tem WR25 was subsequently identified as an SBl by 
Gamen et al. (2006). These statistics undoubtedly 
represent lower limits to the WR binary frequency, 
although one would expect a lower binary frequency 
than for O stars since some WR stars will be merger 
products and others will become single following the 
(dynamical or post-supernova) disruption of initial 
binary systems. 

2 WR stars and star clusters 

If we adopt a lower mass limit to the production of 
WR stars of ^25 Mq, and if massive stars only form 
in clusters, whose upper mass limit follows the rela¬ 
tionship of Weidner et al. (2010), one would expect 
IO^Mq clusters to host the majority of WR stars. 
In reality only 27% of WR stars in the Milky Way 
are in known star clusters (Rosslowe & Growther, 
2015). Of these, the majority of known WR stars in 
the Gentral Molecular Zone (GMZ) do lie within one 
of the three massive star clusters, the Arches, Quin¬ 
tuplet and the Galactic Gentre cluster itself. Ex¬ 
cluding the GMZ, only 18% of the WR stars located 
in the disk lie within star clusters. This is illustrated 
in Fig. 2. 

Gonsequently, either the lower mass limit to the 
production of WR stars must be significantly lower 
(making birth clusters harder to identifity), WR 
stars do not preferentially form in dense star clus¬ 
ters, or the majority of WR stars are no longer as¬ 
sociated with their birth cluster. 

OB stars are observed in a range of environments 


- low density star forming regions (e.g. p Oph), 
intermediate density OB associations such as Gyg 
OB2, and dense clusters such as NGG 3603. Indeed, 
the fraction of stars (of all masses) located in dense 
(Orion Nebula-like) clusters in the Solar Neighbour¬ 
hood is <26% (Bressert et ah, 2010). Wright et al. 
(2014) have established that massive stars in Gyg 
OB2 did not form in close proximity, nor in regions 
of higher density, so can form in regions of relatively 
low density. 

Therefore, it is likely that typical massive stars 
arise from relatively loose OB associations, with only 
a small fraction born in dense star clusters, so the 
rarity of WR stars in star clusters naturally follows 
without resorting to a low threshold mass to the for¬ 
mation of WR stars or a high ejection frequency. 
This has relevance to the discussion of Smith & 
Tombleson (2015) regarding the relative masses of 
O, WR stars and Luminous Blue Variable (LBV) in 
the Large Magellanic Gloud from their spatial loca¬ 
tions. One would expect relatively few O stars in 
close proximity to most WR stars/LBVs if the ma¬ 
jority of massive stars originate in intermediate den¬ 
sity environments, especially since most WR stars 
will arise from relatively modest ^ 25M0 progeni¬ 
tors. 


3 Distribution of WR stars 

Since the distance of the majority of Galactic WR 
stars is not well established, Rosslowe & Growther 
(2015) have investigated the near-IR absolute mag¬ 
nitudes of different WR subtypes based on 108 stars 
whose distances have been estimated, either from 
cluster/association membership or other techniques. 
The spatial distribution of these, plus 246 field WR 
stars for which reliable spectral types have been ob¬ 
tained, are presented in Fig. 3. 

Field WR populations include binary systems, 
providing the WR to OB light ratio can be estimated 
in the K-band or the WR star is thought to domi¬ 
nate the near-IR continuum flux. From a census of 
local (<3 kpc) WR systems, the WR star dominates 
the K-band flux in 82% of cases. In addition, dusty 
WG stars are omitted from this distribution because 
their K-band flux is typically dominated by hot dust 
rather than the WG continuum. For reference, dusty 
WG stars comprise ^15% of the local (<3 kpc) WR 
population. 

Generally, the near-IR classification of WR stars is 
rather coarser than at visible wavelengths. The ma¬ 
jority of near-IR studies have followed the scheme 
of Growther et al. (2006) who utilised IRTF/SpeX 
1-5/im observations of ^30 optically classified WR 
stars. In general robust near-IR classification re¬ 
quires observations at J, H and K, although the ac¬ 
tual number of WN/G transition stars will be un¬ 
derestimated since no near-IR schemes have so far 
been implemented (efforts are currently underway 


2 




WR content of the Milky Way 


to remedy this deficiency). 

# Non-cluster members (excl GC) # Cluster members (excl GC) 



# Non-cluster members (GC) # Cluster members (GC) 



# Non-cluster members # Cluster members 


The observed ratios of WR subtypes in the inner 
Galaxy at galactocentric distances of below 6 kpc 
(log O/H+12 ^ 8.85), mid disk from 6-9 kpc (log 
O/H+12 ^8.7) and outer disk beyond 9 kpc (log 
0/H + 12 ^ 8.55) is presented in Table 2. Overall 
there is little variation between these regions with 
the exception of early to late WC stars, with the lat¬ 
ter dominating in the inner disk owing to the metal- 
licity dependence of WC classification diagnostics 
(Crowther et ah, 2002). Incorporating results from 
the Magellanic Clouds, the observed WC to WN ra¬ 
tio is presented in Fig. 4, together with a variety 
of evolutionary predictions from single and binary 
models (see caption). 


0 -□ 



Large symbols: known distonce 

Small symbols: distance using calibration 



0 

x[kpc 



Fig. 2 : Pie charts illustrating the fraction of Wolf-Rayet 
stars (vl.l3 of catalogue) which are cluster members 
(green) or in the field (blue) for the Galactic disk (top). 
Galactic Centre region (middle) and combined pop (bot¬ 
tom). 


Fig. 3 : Distribution of 354 WR stars in the Galactic 
disk from Rosslowe & Crowther (2015) 


4 Global WR content 

Various estimates of the global WR population in 
the Milky Way have been made, ranging from 1,200 
(Maeder & Lequeux, 1982) to 6,500 (van der Hucht, 
2001). We have constructed a toy model of the 
WR population in the Milky Way in an azimuthally 
symmetric disk following the radial HII distribution, 
atomic/molecular dust distribution and the observed 
WN/WC distribution (Rosslowe & Crowther, 2015). 
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Tab. 2: WR subtype distribution in the Milky Way (Rosslowe & Crowther, 2015) for three galacto-centric distances 
(Rgc)- 


Region 

N(WR) 

N(WC)/N(WN) 

N(WC+d)/N(WN) 

N(WCE)/N(WCL) 

Inner {Rqc < 0 kpc) 

187 

0.51 

0.69 

0.05 

Mid (6 < Rqc < 9 kpc) 

132 

0.53 

0.73 

1.0 

Outer {Rgc > 9 kpc) 

35 

0.40 

0.57 

1.5 


If we assume that the observed WR distribution is 
complete to K=8 mag (though see below), ^550 WR 
dominated stars in the Galactic disk reproduce the 
observed WR distribution. In addition, over 100 WR 
stars are currently known with —2.5° < I < +3.5° 
so the total CMZ population may be as large as 250. 
Such a large population is surprisin given that the 
CMZ accounts for perhaps 5% of the entire Milky 
Way star formation rate (Longmore et ah, 2013), 
but the presence of three very massive clusters in 
the Galactic Centre region suggests WR stars may 
be disproportionately represented. 



Fig. 4: Ratio of WC to WN stars in the Milky Way, 
for which dusty WC stars have been omitted (inverted 
open triangles) or included (inverted filled triangles), 
plus the Magellanic Clouds (open triangles) from Ross¬ 
lowe & Crowther (2015). Predictions from rotating sin¬ 
gle star models (Meynet & Maeder, 2005, green), binary 
models (Eldridge et ah, 2008, red) and non-rotating sin¬ 
gle star models for various mass-loss metallicity depen¬ 
dencies (Eldridge & Vink, 2006, blue) are also indicated. 


A total of 550 + 250 = 800 stars would repre¬ 
sent only 82% of the non-dusty WR population since 
18% of systems are expected to be dominated by the 
WR companion in the K-band, so the non-dusty WR 
population inferred is ^950. Finally, including dusty 
WC stars, which comprise 15% (150) of the whole 
population, we estimate N(WR) = 1,100, providing 


N(WCd)/N(WC) is uniform across all galactocentric 
radii. A histogram of the predicted JHKg magni¬ 
tudes for this predicted WR population is presented 
in Fig. 5. 



Fig. 5: Histogram of JHKs and GAIA G-band magni¬ 
tudes predicted from the model WR population by Ross¬ 
lowe V Crowther (2015).Thick lines represent popula¬ 
tions consisting of WN and non-dusty WC stars, while 
the dotted-line illustrates a model population in which 
28% of WC stars are dust-forming with Mr = —6.95 
mag. 


Based on this global WR population of 1,100, if 
the lower mass threshold to WR stars is 25M0, 
an average lifetime for the WR phase of 0.25 Myr 
is implied by the current Milky Way star forma¬ 
tion rate of 1.9 Mq yr“^ (Chomiuk & Povich, 
2011) and a standard Kroupa Initial Mass Func¬ 
tion. This characteristic lifetime is in fair agre- 
ment with non-rotation evolutionary models, but is 
somewhat shorter than rotating model predictions 
(Georgy et ah, 2012). 

If the observed WR distribution is not complete 
to K=8 mag, the acual WR population will increase. 
By way of example, the two most recent additions to 
the on-line catalogue (vl.l3) involve WRlll-13^, a 
WN6 star from Messineo et al. (2015) with K=8.02 


^All WR stars discovered between the 6th catalogue and the Annex to the 7th catalogue follow the WRXXXa, WRXXXb 
nomenclature, whereas more recent additions utilise WRXXX-1, -2 
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mag and WRlll-14, a WN7-8 star from Nebot 
Gomez-Moran et al. (2015) with K=7.7 mag. If the 
WR-dominated disk population were 650 instead of 
550, a global population of 1,300 would be inferred, 
so the current census may be 50% complete. 

5 Summary 

A brief overview of the Wolf-Rayet population in the 
Milky Way is presented. The current census of 642 
(vl.l3 of online catalogue) represents a doubling of 
the population since the Annex to the 7th WR cata¬ 
logue only a decade ago, due to systematic searches 
from narrow-band surveys at 2/im and broad-band 
near to mid-IR surveys. 

A strict lower limit to the binary fraction is 34% 
(SB 1-2 systems plus dusty WC stars) for V< 15 
mag, but the actual fraction will oudoubedly be sig¬ 
nificantly higher. Spectroscopic monitoring of WR 
stars in the near-IR is now feasible owing to multi¬ 
object IR integral field units (e.g. KMOS at ESO’s 
VLT). 

Excluding the Galactic Centre region, <20% of 
WR stars are located in star clusters, adding weight 
to the formation of most massive stars in OB asso¬ 
ciations, i.e. away from dense clusters. The current 
WR census may be 50% complete, in which case the 
average duration of the WR phase is 0.25 Myr for a 
lower mass limit of 25 Mq to the formation of WR 
stars. 

The majority of systems hosting WR stars are 
dominated by the He-rich component in the K- 
band, so these are fundamentally high luminosity 
(high mass) stars. In contrast, typical progenitors of 
stripped envelope core-collapse supernovae are rela¬ 
tively low mass stars, having transferred the major¬ 
ity of their H-rich envelope to a close companion. 
Such systems would not be recognised as conven¬ 
tional WR stars since the He-rich mass donor would 
be masked by the H-rich mass gainer at optical and 
near-IR wavelengths. 

Einally, the online WR catalogue is maintained in 
Sheffield on a best efforts basis so authors are en¬ 
couraged to provide Paul Crowther with preprints 
as soon as they are accepted for publication. 

References 

Bressert, E., Bastian, N., Gutermuth, R., et al. 2010, 

MNRAS, 409, L54 

Gampbell, W. W. 1894, Astronomy and Astro- 

Physics (formerly The Sidereal Messenger), 13, 
448 

Ghene, A.-N., Borissova, J., Glarke, J. R. A., et al. 

2012, A&A, 545, A54 


Ghomiuk, L. & Povich, M. S. 2011, AJ, 142, 197 
Growther, P. A. 2007, ARA&A, 45, 177 
Growther, P. A., Dessart, L., Hillier, D. J., Abbott, 
J. B., & Eullerton, A. W. 2002, A&A, 392, 653 
Growther, P. A., Hadfield, L. J., Glark, J. S., 
Negueruela, L, & Vacca, W. D. 2006, MNRAS, 
372, 1407 

Eldridge, J. J., Izzard, R. G., & Tout, G. A. 2008, 
MNRAS, 384, 1109 

Eldridge, J. J. & Vink, J. S. 2006, A&A, 452, 295 
Eiger, D. E., Najarro, E., Gilmore, D., et al. 2002, 
ApJ, 581, 258 

Eleming, W. P. S. & Pickering, E. G. 1912, Annals 
of Harvard Gollege Observatory, 56, 165 
Gamen, R., Gosset, E., Morrell, N., et al. 2006, 
A&A, 460, 777 

Georgy, G., Ekstrom, S., Meynet, G., et al. 2012, 
A&A, 542, A29 

Hadfield, L. J., van Dyk, S. D., Morris, P. W., et al. 
2007, MNRAS, 376, 248 

Longmore, S. N., Bally, J., Testi, L., et al. 2013, 
MNRAS, 429, 987 

Maeder, A. & Lequeux, J. 1982, A&A, 114, 409 
Mauerhan, J. G., Van Dyk, S. D., & Morris, P. W. 
2011, AJ, 142, 40 

Messineo, M., Glark, J. S., Eiger, D. E., et al. 2015, 
ApJ, 805, no 

Meynet, G. & Maeder, A. 2005, A&A, 429, 581 
Nebot Gomez-Moran, A., Motch, G., Pineau, E.-X., 
et al. 2015, MNRAS, 452, 884 
Payne, G. H. 1930, Harvard Observatory Mono¬ 
graphs, 3, 1 

Roberts, M. S. 1962, AJ, 67, 79 
Rosslowe, G. K. & Growther, P. A. 2015, MNRAS, 
447, 2322 

Sana, H., de Mink, S. E., de Koter, A., et al. 2012, 
Science, 337, 444 

Shara, M. M., Eaherty, J. K., Zurek, D., et al. 2012, 
AJ, 143, 149 

Shara, M. M., Moffat, A. E. J., Smith, L. E., et al. 
1999, AJ, 118, 390 

Smith, L. E. 1968, MNRAS, 138, 109 
Smith, N. & Tombleson, R. 2015, MNRAS, 447, 598 
van der Hucht, K. A. 2001, New A Rev., 45, 135 
van der Hucht, K. A. 2006, A&A, 458, 453 
van der Hucht, K. A., Gonti, P. S., Lundstrom, L, & 
Stenholm, B. 1981, Space Sci. Rev., 28, 227 
Wachter, S., Mauerhan, J. G., Van Dyk, S. D., et al. 
2010, AJ, 139, 2330 

Weidner, G., Kroupa, R, & Bonnell, I. A. D. 2010, 
MNRAS, 401, 275 

Wright, N. J., Parker, R. J., Goodwin, S. P., & 
Drake, J. J. 2014, MNRAS, 438, 639 


5 




